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ABSTRACT 

The  direction  of  flow  along  the  centerline  of  several  sections 
of  the  inlet  duct  of  a  centrifugal  pump  was  determined  using  a  hot  wire 
anemometer.  In  contrast  to  other  tests  of  this  kindj,  the  pump  was 
operated  with  air  to  facilitate  the  measurements .  The  present  investi- 
gation is  part  of  a  program  at  the  U.  S.  Naval  Postgraduate  School  to 
evaluate  the  performance  of  liquid  pumps  with  air  in  order  to  apply 
aerodynamic  measuring  methods  and  equipment  to  the  task„  Data  were 
taken  with  the  pump  operating  at  different  points  on  the  characteristic 
curve.  The  flow  direction  data  are  presented  as  a  function  of  position 
in  the  intake  duct  for  various  operating  conditions. 

Prerotation  was  observed,  the  character  of  which  was  in  agreement 
with  theoretical  predictions  in  that  prerotation  was  in  the  same  direc- 
tion as  the  impeller  rotation  at  the  low  flow  rates  and  in  the  opposite 
direction  at  the  high  flow  rates.  Backflow  near  the  periphery  of  the 
duct  was  observed  at  the  low  flow  rates. 
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AN  INVESTIGATION  OF  THE  PREROTATION  CHARACTERISTICS 
OF  A  FLUID  IN  THE  INTAKE  DUCT  OF  A  CENTRIFUGAL  PUMP 

1.  Introduction . 

It  has  been  observed  for  many  years  that  rotation  of  fluid  in  the 
inlet  ducts  ahead  of  an  impeller  exists .  This  rotation  may  be  in  the 
same  direction  as,  or  opposite  to,  the  direction  of  rotation  of  the 
impeller.  The  effect  is  more  pronounced  and  easily  observable  at  off 
design  conditions  but  may  even  occur  at  the  design  point 0  It  is  essen- 
tial to  know  the  fluid  angles  at  the  impeller  inlet  in  order  to  properly 
design  the  inlet  blade  angles.  If  the  blade  angle  is  very  different 
from  the  relative  flovr  angles,  local  high  fluid  velocities  may  produce 
the  phenomenon  of  cavitation „  While  some  cavitation  may  be  tolerated 
to  achieve  higher  flow  rates,  it  is  generally  undesirable  because  it 
is  destructive  and  causes  inefficient  and  unstable  operation  of  the 
pump„ 

The  cause  of  the  phenomenon  o  prerotation  is  somewhat  obscure  and 
little  quantitative  data  is  available  to  pump  designers „  Although  con- 
ventional hydraulic  pumps  have  been  examined  for  prerotation,  there  do 
not  exist  prerotation  data  in  the  literature  for  the  very  small  blade 
angles  encountered  in  rocket  fuel  and  oxidizer  pumps.  Experiments  that 
have  been  done  generally  have  used  liquid  pumps  with  water  as  the  working 
fluid.  It  will  be  seen  that  the  same  results  may  be  obtained  using  a 
liquid  pump  but  with  air  as  the  working  fluid.  The  latter  would  require 
much  less  power  to  drive  it,  less  elaborate  experimental  equipment,  and 
would  provide  easier  ways  to  investigate  the  physical  aspects  of  the  process, 

2.  The  Pump. 

The  pump  used  for  this  experiment  was  a  liquid  fuel  pump  for  the 
Titan  I  missile.  It  is  a  centrifugal  pump  seven  inches  in  diameter  at 


the  eye  and  11  inches  at  the  discharge  of  the  impeller,,  It  has  a  diffuser 
with  four  vanes,  an  annular  collector  and  a  three  inch  diameter  discharge 
flange.  There  is  a  static  pressure  tap  at  the  discharge  flange  where 
the  pressure,  ?2»   was  "taken .  The  impeller  has  two  blades  at  the  eye 
with  additional  blades  starting  at  larger  radii  for  a  total  of  eight 
blades  at  the  discharge.  Figs.  1,  2,  3.  The  pump  was  driven  by  an  air 
turbine . 
3.  The  Inlet  Duct. 

The  inlet  duct  consisted  of  a  seven  inch  diameter ,  four  inch  long 
plexiglas  spacer,  a  seven  inch  diameter,  2U   inch  long  steel  tube  and 
a  built  up  wooden  bell  mouth.  Through  the  plexiglas s  qualitative 
observations  of  the  flow  were  made  using  smoke  and  cotton  tufts.  Mounted 
on  the  steel  tube  were  four  brackets  with  a  hole  through  the  pipe  'kallj 
through  which  a  probe  could  be  inserted  to  make  a  survey  of  the  flow 
along  the  centerline  of  the  pipe,  top  to  bottom.  These  survey  points 
were  six  inches  apart,  with  the  nearest  one  to  the  impeller  being  seven 
inches  from  the  inlet  flange.  This  point  is  referred  to  as  position 
number  four. 
4-.  The  Discharge  Pipe. 

The  discharge  pipe  was  a  three  inch  pipe,  ten  feet  long.  It  had 
a  bundle  of  one -half  inch  copper  tubing  six  inches  long  at  the  pump  end 
to  act  as  a  flow  straightener.  Five  feet  from  the  pump  was  a  flow  meter 
consisting  of  a  total  pressure  Kiel  probe  with  a  thermocouple  mounted  at 
the  centerline  of  the  pipe  and  two  static  ports  on  opposite  sides  of  the 
pipe.  For  calibration  of  the  flow  meter  see  Appendix  C.  At  the  discharge 
end  of  the  pipe  was  a  gate  valve  for  controlling  the  flow  rate. 


5.  The  Hot  Wire  Anemometer. 

A  Flow  Corporation,  Model  HWB,  hot  wire  anemometer  was  used  to 
determine  flow  direction  using  a  O00035  inch  diameter  tungsten  wire. 
The  probe  was  mounted  in  a  protractor  which  in  turn  was  mounted  on  a 
slide.  The  depth  of  the  hot  wire  could  be  determined  by  a  scale  on 
the  slide  and  the  directional  orientation  of  the  hot  wire  could  be  read 
to  0.2  degrees  using  a  vernier  scale  on  the  protractor.  The  centerline 
of  the  inlet  duct  at  position  number  two  at  the  highest  flow  rate  was 
used  as  a  calibration  point  for  directionally  orienting  the  hot  wire. 

6.  Other  Equipment. 

An  electronic  counter  was  used  to  determine  pump  speed.  A  bank  of 
vertical  tube  manometers  was  used  to  measure  pressures  in  cm.  of  water. 
The  differential  pressure  was  measured  with  a  water  micromanometer0 
Discharge  total  temperature  was  determined  with  a  Brown  thermocouple 
instrument.  A  four -foot  cubic  plywood  box  was  used  as  an  inlet  plenum 
to  act  as  a  buffer  between  the  inlet  duct  and  ambient  turbulence.  The 
inlet  to  the  plenum  was  a  two-foot  square  hole  covered  with  four  layers 
of  filter  cloth. 

7.  Speed  of  the  Pump. 

The  design  speed  of  the  pump  is  14-, 000  RPM.  It  was  originally 
planned  that  several  speeds  in  the  neighborhood  of  the  design  speed 
would  be  used  for  this  experiment.  Upon  running  the  pump  at  speeds  of 
10,000  to  12,000  RPM,  large  amplitude  vibrations  were  observed.  A 
dynamic  balance  of  the  air  turbine  failed  to  lessen  the  vibrations  so 
a  speed  of  8,800  RPM  was  settled  upon  for  the  tests.  The  vibrations 
are  evidently  due  to  the  ball  bearings  whose  axial  loading  is  not  as 
great  as  design  loading  because  air  is  the  fluid  being  pumped. 


In  order  to  avoid  cavitation  the  axial  velocity  ahead  of  the  pump 
is  kept  small.  The  reduction  in  RPM  produced  a  further  decrease  in 
velocity  which  caused  difficulty  in  measuring  the  dynamic  pressure  in 
the  duct, 
8.  Measuring  the  Flow  Angle. 

The  first  attempt  at  measuring  the  flow  direction  in  the  inlet  duct 
was  made  utilizing  a  United  Sensor  airflow  probe.  This  probe  had  five 
pressure  ports.  When  the  probe  was  properly  oriented  in  the  air  stream, 
one  would  measure  the  total  pressure,  two  would  measure  the  static 
pressure,  and  two  others  not  used  in  this  application  could  be  used  to 
measure  components  of  the  flow  in  the  direction  of  the  probe  axis. 
After  calibration  in  an  air  stream  of  known  direction,  the  probe  could 
be  used  to  determine  the  direction  of  airflow  in  the  inlet  duct.  Not 
only  the  direction,  but  also  the  velocity  of  the  flow  could  be  determined 
at  any  point  in  the  duct  tested,  and  from  this  the  flow  rate  could  have 
been  calculated. 

The  very  low  flow  rates  made  the  pressure  measurements  to  determine 
flow  angles  very  critical.  The  maximum  flow  rate  which  was  measured 
later  with  a  flow  meter  was  .4.6  lb/sec.  This  flow  rate,  assuming  a 
uniform  velocity  distribution  in  the  inlet  duct,  would  yield  a  dynamic 
head  of  about  0.1  inches  of  water.  Consider  now  the  ideal  flow  about 
a  cylinder  where 

Pg  -  Pt  :  q0(l-Asin2e)  =  0.l(l-4sin2e) 
then   dPg  B  -0.8  sin@cosQ 

de 

for     6   =  45° 

dPg_  =  -0.4.  in.H20/radian  -  -0.007  inoH20/degree 

d0 


Hence,  if  the  probe  were  one  degree  away  from  facing  directly  into  the 
stream,  the  pressure  difference  between  the  two  static  ports  would  be 
about  „015  inches  of  water  <>  An  inclined  manometer  was  available  which 
could  be  read  to  ,01  inches „  However,  since  the  sensitivity  goes  down 
even  further  as  the  flow  rate  is  reduced,  it  was  felt  that  the  accuracy 
of  results  would  be  rather  poor  with  this  method 0 

The  hot  wire  anemometer  offered  a  possible  solution „  As  pointed 
out  in  Appendix  D,  the  directional  sensitivity  of  the  hot  wire  is  great- 
est when  the  wire  axis  is  aligned  with  the  direction  of  flowu  It  was 
felt  that  the  supporting  pin  ahead  of  the  wire  might  cause  turbulence 
and  shield  the  flow  from  the  full  effect  of  the  cooling  gas  stream „  To 
avoid  this,  the  leading  end  support  pin  was  filed  off  so  that  the  wire, 
when  soldered  across  the  tips  of  the  pins,  made  an  angle  of  about  65° 
with  the  probe  axis0  No  apparent  reduction  in  the  galvanometer  needle 
fluctuations  was  noted ,  Since  the  directional  sensitivity  of  the  probe 
is  reduced  by  this  technique,  the  standard  probe  was  used. 

The  hot  wire  anemometer,  being  sensitive  enough  to  detect  very 
small  velocities,  is  also  sensitive  to  perturbation  velocities „  There- 
fore, it  became  necessary  to  build  the  intake  plenum  chamber,  described 
in  Sec0  6,  to  prevent  ambient  turbulence  from  giving  an  indication  of  a 
change  in  flow  direction .  Even  with  the  plenum  chamber  there  was  always 
some  fluctuation  of  the  galvanometer  needle  used  to  determine  the  resis- 
tances „  With  fluctuations  of  the  needle  it  was  difficult  to  find  a 
maximum,  so  the  scheme  actually  used  was  to  record  the  two  angles, 
15-20  degrees  on  either  side  of  the  angle  for  maximum  deflection,  which 
produced  the  same  galvanometer  reading .  The  angle  halfway  between  the 
two  recorded  was  then  used  as  the  flow  angle „  In  this  technique  there 


can  be  no  shielding  by  the  upstream  supporting  pin  because  it  is  not 
directly  in  front  of  the  wire0 
9.  Existing  Theories „ 

While  the  exact  mechanism  of  prerotation  is  not  clear,  the  phenom- 
enon has  been  observed  by  many  investigators  for  years  .  0>J°  Several 
attempts  have  been  made  to  explain  the  reasons  for  prerotation „  At  c$fcf 
design  conditions  the  fluid  is  said  to  acquire  prerotation  by  taking  the 
path  of  least  resistance  or  least  inlet  impact.,  For  flow  rates  lower 
than  design  the  flow  acquires  prerotation  which  is  the  same  direction 
as  the  impeller  rotation,  a  positive  prerotation,  and  at  flow  rates 
higher  than  design  the  prerotation  is  negative  0  jjff »  Backflow  in  the 
outer  portion  of  the  impeller  channel  at  low  flow  rates,  and  high  net 
heads  are  believed  to  be  the  cause  of  positive  prerotation0  Backflow 
has  been  observed  in  the  inlet  duct  and  in  the  diffuser 0  {Cj  „ 

The  impeller  discharge  conditions  seem  to  have  a  great  deal  of 
influence  on  prerotation .  Experiments  at  California  Institute  of 
Technology  using  a  vaneless,  symmetrical  diffuser  showed  practically 
no  evidence  of  prerotation  over  the  entire  range  of  flow  rates „ [2j  0 
This  pump  did  show  a  positive  component  of  tangential  velocity  at  the 
discharge  of  the  impeller  when  operating  at  zero  pressure  coefficient., 
This  positive  angular  momentum  at  the  discharge  would  seem  to  indicate 
a  positive  rotation  at  the  inlet,  but  no  such  rotation  was  observable „ 
At  zero  pressure  coefficient  the  accepted  theory  would  predict  a  negative 
prerotation.  Shepherd  believes  that  the  positive  angular  momentum  at  the 
discharge  is  due  to  losses  in  the  impeller.,  t5j° 

The  pump  studied  by  Peck  had  a  vaneless  diffuser,  but  had  the  typi- 
cal eccentric  discharge  volute  „  j[4]  u  This  pump  exhibited  considerable 


prerotation  and  backflow  at  low  flow  rates „  It  showed  no  negative 
prerotation  at  flow  rates  above  the  design  point e 

Existing  theories,  therefore ,  seem  to  be  of  little  practical  use 
to  the  pump  designer0  The  discharge  conditions  evidently  are  a  contrib- 
uting influence  on  prerotation „ 
10.  Results 8 

The  results  of  these  tests  are  tabulated  in  Tables  1  through  22 
and  are  presented  in  graphical  form  in  Figs0  8  through  17 „ 
11 o  Discussion ■ 

Fig.  8  is  the  experimentally  determined  characteristic  curve  of 
the  pump0  The  pressure  coefficient  used  here  is  non-dimensional  and 
is  defined  as  follows  % 

sy2   =AP//52U22 
where,  p2   is  a  derived  density  using  the  static  pressure  at  the  outlet 
flange  and  the  total  temperature  at  the  flow  meter u  The  flow  coefficient 
is  also  non-dimensional  and  is  defined  as  follows i 

♦2  =  W&2 

Vm2,  not  being  measured  directly,  is  a  derived  quantity  using  the 

measured  flow  rate,  w,  the  density,  p%%   and  outlet  area  of  the  impeller 0 
The  outlet  area  does  not  consider  any  restriction  by  the  blades  and  is 
merely  TfD2b0 

Vm2  =  ^2TnrD2b 
The  theoretical  pressure  coefficients  as  a  function  of  the  flow 

coefficient  were  computed  using  the  slip  factors  of  Stodola  and  Busemann0 

These  are  plotted  in  Fig„  6  for  reference «  The  Stodola  slip  factor^ jaa9 

was  calculated  using  his  equations  \[5) 

pa  s    1  -  IT  cos^2B       =     \2 

n  Vu2  (theoretical) 


Busemann's  slip  factor, yug,  for  eight  blades  and  an  exit  blade  angle 

of  -73.750>  Fig.  3,  was  obtained  from  a  graph  in  Ref0  5.  The  theoretical 

pressure  coefficient  then  may  be  expressed  ass 

The  efficiency  may  be  expressed  as  the  ratio  of  actual  to  theoret- 
ical pressure  coef f icients „  The  efficiency  was  calculated  and  plotted 
in  Fig.  7.  The  pressure  coefficients  and  efficiency  can  be  expected  to 
increase  somewhat  with  the  liquid  fuel  as  a  working  fluid  due  to  increased 
Reynolds  number., 

The  prerotation  observed  in  these  tests  is  generally  what  would  be 
predicted  by  accepted  theories.  Positive  prerotation  and  baekflow  in 
the  outer  portions  of  the  inlet  duct  were  observed  at  low  flow  rates 0 
Positive  prerotation  is  indicated  in  the  figures  as  an  angle  greater 
than  zero  for  depths  less  than  three  and  one  half  inches  and  an  angle 
less  than  360  degrees  for  greater  depths 0 

There  is  considerable  scatter  in  the  data  with  respect  to  the  flow 
angle.  Several  repeatability  tests  were  made  during  the  course  of  the 
experiment  which  indicated  that  the  data  should  be  accurate  to  within 
one  degree.  These  tests  were  performed  in  connection  with  the  direc- 
tional calibration  of  the  hot  wire.  Each  time  the  wire  had  to  be 
replaced  it  was  necessary  to  establish  a  new  angular  tare.  Since  this 
was  a  one  time  measurement  which  effected  all  the  readings  taken  with 
that  wire,  it  was  considered  worthwhile  to  use  an  average  of  eight  to 
ten  readings  rather  than  just  one.  No  readings  were  more  than  one  degree 
different  from  the  average.  However,  these  tests  were  made  on  the 
centerline  of  the  inlet  duct  at  position  number  two  with  the  pump 
operating  at  maximum  flow  coefficient.  Directional  sensitivity  of  the 
hot  wire  is  proportional  to  the  velocity  so  that  the  accuracy  would  be 
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poorer  at  the  lower  flow  rates .  Velocity  perturbations  in  any  direction 
would  give  an  indication  of  a  change  in  flow  direction  to  the  observer. 
It  can  reasonably  be  expected  that  these  perturbations  would  be 
oscillatory  in  character,  however , 

Increasing  the  speed  of  the  pump  to  14. ,000  RPM  would  considerably 
increase  the  accuracy  of  the  data.  Directional  sensitivity  is  propor- 
tional to  the  velocity  with  the  hot  wire  anemometer  and  proportional 
to  the  square  of  the  velocity  if  the  pi tot-static  tube  is  used.  Velocity 
is  proportional  to  the  RPM.  Perturbation  velocities  would  probably  not 
increase  as  much  so  the  fluctuations  of  the  needle  would  be  reduced  also. 

In  spite  of  the  scatter,  it  is  felt  that  some  negative  prerotation 
was  present  at  the  high  flow  rates.  In  Fig.  17,  an  interpretation  has 
been  made  which  shows  a  maximum  of  two  degrees  of  negative  prerotation 
at  the  maximum  flow  rate. 

At  the  apparent  design  point  there  is  still  considerable  positive 
prerotation  and  backflow.  The  point  of  highest  efficiency  is  at  a  flow 
coefficient  of  .03.  Fig.  12,  showing  data  taken  at  a  flow  coefficient 
of  .04.8,  shows  that  positive  prerotation  and  backflow  were  observed  at 
position  number  four. 

Velocity  measurements  were  not  taken  so  the  exact  incidence  angles 
cannot  be  calculated.  However,  by  assuming  that  the  magnitude  of  the 
velocity  is  a  constant  across  the  section  for  each  flow  rate,  an  estimate 
of  the  effect  of  prerotation  on  the  incidence  angles  may  be  made.  This 
was  done  for  the  two  extreme  flow  rates.  The  following  analysis  was 
made  to  determine  the  velocities  in  the  inlet  duct  for  the  case  with 
prerotation  and  for  the  case  of  axial  flow  with  uniform  velocity 
distribution  for  comparison. 


Vs*  -  ^2U2A2  =  2TTv/r  cos  c<x  dr 

for  the  axial  flow  case 
(f)2U2A2  =  T?fT(r02  -  r±2) 
The  results  weres 

t2   =  .026      V(wlth  rotation)  =  33.6  ft/sec 

V(  axial  flow)    =   7.1  ft/sec 

*2  =  *097      V(with  rotation)  =  31»2  ft/*ec 

v( axial  flow)    -  26.-4  ft/sec 

These  results  were  used  in  plotting  the  velocity  triangles  in  Figs.  19 

through  22 . 

It  can  be  seen  that,  except  for  the  very  large  amounts  of  prero- 
tation  at  the  low  flow  rate,  velocity  of  the  flow  has  a  greater  influence 
on  the  incidence  angles  than  the  flow  angle  itself „  Of  course,  these 
velocities  must  be  greater  with  prerotation  than  without  for  any  given 
flow  rate.  Therefore,  the  tendency  to  cavitate  would  be  greater  with 
prerotation  because  the  increased  velocities  produce  reduced  static 
pressures. 
12 ,  Conclusions. 

These  tests  show  essentially  the  same  trands  as  other  tests  have 
shown  where  water  was  used  as  the  working  fluid. 

The  very  small  velocities  involved  may  be  handled  very  well  with 
the  hot  wire  anemometer.  This  instrument  may  be  used  to  measure  the 
velocity  as  well  as  the  flow  direction.  Appendix  D0 

Testing  liquid  pumps  with  air  has  many  advantages.  The  power 
requirements  for  driving  the  pump  are  less.  Return  and  supply  lines 
with  their  accompanying  pumps  and  coolers  are  not  required.  Measuring 
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equipment  is  not  subjected  to  high  dynamic  loads „  Test  setups  may  be 
kept  simple  and  aerodynamic  measuring  devices  can  be  used*,  Packing  is 
not  required  around  the  probe  to  prevent  leakage  and  the  unused  holes 
in  the  inlet  duct  may  be  sealed  with  a  piece  of  drafting  tape0 
13.  Recommendations . 

It  is  recommended  thats 

a„  Efforts  be  made  to  eliminate  vibrations  so  that  tests 
may  be  made  at  higher  speeds  0  Preloading  the  bearings 
and  a  dynamic  balance  of  the  entire  assembly  should  reduce 
vibrations o  Installation  in  one  of  the  new  test  cells  in 
a  more  permanent  fashion  would  eliminate  feedback  of 
vibrations  from  the  supply  line  and  the  discharge  pipec 
Directional  sensitivity  of  the  hot  wire  is  proportional 
to  the  velocity  so  this  would  improve  the  accuracy  of  the 
data„ 
b0  Velocity  measurements  as  well  as  direction  measurements 
be  taken  using  the  hot  wire  anemometer*  Direction  should 
be  determined  first  and  then  velocity  should  be  measured 
with  the  wire  axis  perpendicular  to  the  air  stream, 
c.  A  more  sophisticated  probe  which  employs  two  hot  wires 
on  the  same  probe  be  purchased  or  fabricated  to  try  to 
obtain  more  accurate  data* 
d„  The  diffuser  vanes  be  removed  to  determine  the  effect  of 
that  section  of  the  pump  on  the  prerotation0 
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APPENDIX 
TABLES 


H 


TV 


TABLE   1 


PREROTATION  TEST 


RUN  NO.   1 


12  MARCH  1964 


POSITION  NO.  1 


FLOW  COEFFICIENT  =  .07594     PRESSURE  COEFFICIENT  =  .27850 


DEPTH 


ANGLE 


.10 

355.40 

.50 

2.60 

1.00 

1  .60 

1.50 

357.10 

2.00 

357.80 

2.50 

358.70 

3.00' 

359.20 

3.50 

.20 

u.oo 

359.60 

4.50 

1  .40 

5.00 

1  .80 

5.50 

.30 

6.00 

3.80 

6.50 

1  .30 

6.90 

356.50 

IS 


T  AR  L  E   2 


PREROTATION  TEST 


RUN  NO.   2 


13  MARCH  1964 


POSITION  NO.  1 


FLOW  COEFFICIENT  =  .09598     PRESSURE  COEFFICIENT  =  .08463 


CEPTH 


ANGLE 


.10 

.90 

.50 

1.20 

1  .00 

2.00 

1.50 

355.60 

2.00 

357.00 

2.50 

359.60 

3.00 

2.70 

3.50 

.00 

4.00 

3.70 

4.50 

359.40 

5.00 

3.80 

5.50 

1  .10 

6.00 

5.50 

6.50 

355.40 

6.90 

355.70 

16 


u 


'   TABLE   3 
PREROTATION  TEST 
RUN  NO.   3  13  MARCH  196U 

POSITION  NO.  U 
FLOW  COEFFICIENT  =  .09590     PRESSURE  COEFFICIENT  =  .08522 

OEPTH  ANGLE 


.10 

1  .90 

.50 

359.00 

1.00 

359.10 

1.50 

353.95 

2.00 

358.50 

2.50 

357.30 

3.00 

359.10 

3.50 

.95 

U.00 

.50 

4.50 

2.15 

5.00 

1.50 

5.50 

1.10 

6.00 

.20 

6.50 

358.70 

6.90 

359.20 

17 


T  A8  L  E       U 


PRESTATION    TEST 


RUN    NO.       4 


13    MARCH    1964 


POSITION    NO.    U 


FLOW  COEFFICIENT  =  .06336 


PRESSURE  COEFFICIENT  =  .37310 


DEPTH 


ANGLE 


.10 

1.60 

.50 

359.70 

1.00 

.45 

1.50 

352.80 

2.00 

357.00 

2.50 

357.20 

3.00 

358.95 

3.50 

1  .50 

U.00 

.35 

4.50 

359.40 

5.00 

358.75 

5.50 

.30 

6.00 

.15 

6.50 

359.90 

6.90 

356.55 

18 


T  AB  L  E   5 


PREROUTION  TEST 


RUN  NO.   5 


13  MARCH  196U 


POSITION  NO.  4 


FLOW  COEFFICIENT  =  .021U5     PRESSURE  COEFFICIENT  =  .56532 


DEPTH 


ANGLE 


.10 

102.50 

.50 

97.00 

1.00 

84.60 

1.50 

67.20 

2.00 

45.90 

2.50 

27.95 

3.00 

12.20 

3.50 

1.30 

u.oo 

353.30 

4.50 

340. U5 

5.00 

322. UO 

5.50 

299. UO 

6.00 

280.20 

6.50 

26U.95 

6.90 

256.10 

19 


T  Aft  L  E   6 


PREROTATION  TEST 


RUN  NO.   6 


17  MARCH  196U 


POSITION  NO.  U 


FLOW  COEFFICIENT  =  .04808     PRESSURE  COEFFICIENT  =  .49974 


CEPTH 


ANGLE 


.10 

101.45 

.50 

95.00 

1.00  ' 

78.25 

1.50 

55.70 

2.00 

21.50 

2.50 

9.50 

3.00 

2.80 

3.50 

360f00 

4.00 

358.00 

4.50 

354.10 

5.00 

341.40 

5.50 

321.00 

6.00 

289.70 

6.50 

267.75 

6.90 

252.30 

20 


N 


TABLE   7 


PREROTATION  TEST 


RUN  NO-   7 


17  MARCH  1964 


POSITION  NO.  4 


FLOW  COEFFICIENT  =  .06138     PRESSURE  COEFFICIENT  =  .140929 


CEPTH 


ANGLE 


.10 

357.50 

.50 

357.60 

1.00 

354.90 

1.50 

355.60 

2.00 

357.30 

2.50 

357.30 

3.00 

356.50 

3.50 

.05 

4.00 

359.65 

4.50 

354.00 

5.00 

359.35 

5.50 

358. UO 

6.00 

359.70 

6.50 

35  5.40 

6.90 

351 .90 

21 


TABLE   8 


PREROTATION  TEST 


RUN  NO.   8 


18  MARCH  1964 


POSITION  NO.  U 


FLOW  COEFFICIENT  =  .07989     PRESSURE  COEFFICIENT  =  .27829 


CEPTH 


ANGLE 


.10 

.15 

.50 

.40 

1.00 

359.50 

1.50 

.10 

2.00 

358.85 

2.50 

359.90 

3.00 

.45 

5.50 

.40 

4.0C 

.65 

4.50 

358.85 

5.00 

359.50 

5.50 

359.95 

6.00 

1.15 

6.50 

2.05 

6.90 

1.55 

22 


?s 


TABLE       9 


PRESTATION    TEST 


RUN    NO.       9  18    MARCH    1964 


POSITION    NO.    4 


FLOW  COEFFICIENT  *  .06496     PRESSURE  COEFFICIENT  =  .40275 


CEPTH  ANGLE 


;io 

3.45 

.50 

2.40 

1.00 

.90 

1.50 

353.80 

2.00 

359.60 

2.50 

359.35 

3.00 

358.90 

3.50 

359.70 

4.00 

1.00, 

4.50 

.50 

5.00 

358.60 

5.50 

1.145 

6.00 

.75 

6.50 

2.80 

6.90 

.90 

23 


TABLE  10 
PRER0TATI0N  TEST 
RUN  NO.  10  •  18  MARCH  196U  , 

POSITION  NO.  3 
FLOW  COEFFICIENT  =  .065U3     PRESSURE  COEFFICIENT  =  .U0C59 

DEPTH  ANGLE 


.10 

2.60 

.50 

3.20 

1.00 

.70 

1.50 

358.55 

2.00 

359. UO 

2.50 

357.05 

3.00 

359.15 

3.50 

359.55 

U.OO 

359.10 

U.50 

1.95 

5.00 

.80 

5.50 

.50 

6.00 

359.60 

6.50 

1  .05 

6.90 

.  .50 

2U 


TABLE     11 

'I 

•  ti 

PREROTATION    TEST 

»' 

RUN    NO.     11  18    MARCH    1964. 

POSITION    NO.     3 
FLOW   COEFFICIENT    =    .09707*  PRESSURE   COEFFICIENT    =    .08362 

DEPTH  ANGLE 


.10 

1.10 

.50 

2.60 

1.00 

.50 

1.50 

358.25 

2.00 

359.10 

2.50 

359.00 

3.00 

359.60 

3.50 

.20 

4.00 

1  .60 

4.50 

1.30 

5.00 

1.35 

5.50 

"359.30 

6.00 

.45 

6.50 

359.30 

6.90 

.75 

25^ 


■I  ». 


TABLE  12 


PRESTATION  TEST 


RUN  NO.  12 


19  MARCH  1964 


POSITION  NO,  1 


FLOW  COEFFICIENT  =  .08693 


PRESSURE  COEFFICIENT  =  .20918 


DEPTH 


ANGLE 


.10 

14.10 

.50 

1  .00 

1.00 

.50 

1.50 

357. 40 

2.00 

2.10 

2.50 

357.00 

3.00 

358.20 

3.50 

359.50 

4. CO 

359.50 

4.50 

.20 

5.00 

.10 

5.50 

359.50 

6.00 

358.10 

6.50 

359.05 

6.90 

358.50 

26 


TABLE  13 


PREROTATION  TEST 


RUN  NO.  13 


24  MARCH  1964 


POSITION  NO.  1 


FLOW  COEFFICIENT  =  .03853 


PRESSURE  COEFFICIENT  =  .56494 


OEPTH 


ANGLE 


.10 

5.70 

.50 

5.50 

1  .00 

3.00 

1.50 

2.20 

2.00 

1.20 

2.50 

355.80 

">' 

3.00 

352.  pO 

3.50 

353.30 

4.  00 

357.80 

4.50 

357.90 

'  5.00 

359.80 

5.50 

.50 

6.00 

.70 

6.50 

359.60 

6.90 

355.30 

27 


T  AB  L  E    1*4 
PREROTATION    TEST 

J 

RUN  NO.  14  31  MARCH  1964 

POSITION  NO.  2, 
FLOW  COEFFICIENT  =  .09675     PRESSURE  COEFFICIENT  =  .08516 

DEPTH  ANGLE 


'  .10 

2.40 

.50 

4.50 

1.00 

2.20 

1.50 

358.00 

2.00 

356.30 

2.50 

358.60 

3.00 

359^70 

3.50 

.00 

4.00 

1.35 

4.50 

»'.20 

5.00 

3.00 

5.50 

1.80 

6.00 

359.40 

6.50 

358.70 

6.90 

357.00 

28 


TABLE    15 


PRESTATION    TEST 


RUN    NO.     15      . 


31    MARCH    1964 


POSITION    NO.    2 


FLOW  COEFFICIENT  =  .07927     PRESSURE  COEFFICIENT  =  .28404 


DEPTH 


ANGLE 


.10 

558.70 

.50 

358.90 

1.00 

357.70 

1.50 

.15 

2.00 

358.  45 

2.50 

358.60 

3.00 

.10 

3.50 

358.60 

4.00 

.80 

4.50 

1  .50 

•  5.00 

1.50 

5.50 

1.30 

6.00 

359.20 

6.50 

358.70' 

6.90 

.45 

29 


TABLE  16 
PREROTMION  TEST 
RUN  NO.  16  1  APRIL  1964 

POSITION  NO.  2 
FLOW  COEFFICIENT  =  .06435     PRESSURE  COEFFICIENT  =  .40481 
■;'•■.  :■'"   DEPTH  ANGLE 


.10 

3.00 

.50 

1  .90 

1.00 

359.35 

1.50 

358.30 

2.00 

358. UO 

2.50 

358.90 

3. CO 

.70 

3.50 

1.05 

4.00 

.40 

4.50 

.75 

5.00 

360.00 

5.50 

1  .50 

6.00 

359.50 

6.50 

.90 

6.90 

1  .30 

30 


TABLE  17 


PRESTATION  TEST 


RUN  NO.  17 


1  APRIL  196U 


POSITION  NO,  2 


FLOW  COEFFICIENT  =  .04756     PRESSURE  COEFFICIENT  =  .52388 


CEPTH 


ANGLE 


.10 

U.90 

.50 

359.60 

1.00 

.60 

1.50 

358.10 

2.00 

357.70 

2.50 

357.90 

3.00 

360.00 

3.50 

359.80 

4.00 

.30 

4.50 

.  40 

5.00 

1.35 

5.50 

.30 

6.00 

.70 

6.50 

359.40 

6.90 

3.30 

31 


TABLE  18 


PREROTATION  TEST 


RUN  NO,  18 


1  APRIL  1964  , 


POSITION  NO.  2 


FLOW  COEFFICIENT  =  .02765     PRESSURE  COEFFICIENT  =  .60983 


DEPTH 


ANGLE 


.10 

10U.30 

.50 

100.30 

1.00 

87.10 

1.50 

68.20 

2.00 

43.10 

2.50 

22.60 

3.00 

1  1.50 

3.50 

2.20 

U..  00 

357.40 

U.50 

346.50 

5.00 

330.50 

5.50 

306.30 

6.00 

282.60 

6.50 

26U. 10 

6.90 

256.60 

32 


T  AB  L  E  19 


PRESTATION  TEST 


RUN  NO.  19 


1  APRIL  1964 


POSITION  NO.  3 


FLOW  COEFFICIENT  =  .02663     PRESSURE  COEFFICIENT  =  .6C978 


DEPTH 


ANGLE 


.10 

103.30 

.50 

96. ao 

1.00 

87.30 

1  .50 

74.20 

2.00 

56.20 

2.50 

i   35.90 

3.00 

18.10 

3.50 

0.70 

4.00 

352.60 

4.50 

334.90 

5.00 

3J.I*60 

5.50 

295.85 

6.00 

277.60 

6.50 

263. MO 

6.90 

256.60 

33 


TABLE    20 


PREROTATION    TEST 


.    i,  * 


RUN    NO.    20 


1    APRIL    1964 


POSITION    NO.    3 

9  • 

FLOW  COEFFICIENT  =  .07792     PRESSURE  COEFFICIENT  =  .29394 


DEPTH 


ANGLE 


.10 

2.40 

.50 

2.10 

1.00 

359.50 

1.50 

359.00 

2.00 

'   358.35 

2.50 

.20 

3.00 

.60 

3.50 

.30 

4.00 

1.30 

4.50 

.60 

5.00 

359.70 

5.50 

.80 

6.00 

359.80 

6.50 

360.00 

6.90 

359.60 

3U 


TABLE  21 


PREROTATION  TEST 


RUN  NO.  21 


2  APRIL  1964 


POSITION  NO.  3 


FLOW  COEFFICIENT  =  .0869U 


PRESSURE  COEFFICIENT  =  .19722 


DEPTH  

ANGLE 

.10 

.10 

.50 

359.80 

1.00 

357.50 

1.50 

358.20 

2.00 

359.00 

2.50 

.60 

3.00 

360.00 

3.50 

.70 

4.00 

1  .00 

4.50 

1.60 

5.00 

1.10 

5.50 

36G.00 

6.0C 

359.50 

6.50 

359.70 

6.90 

.20 

3* 


TA3LE  22 
PREROTATION  TEST 
RUN  NO.  22  2  APRIL  1964 

POSITION  NO.  4 
FLOW  COEFFICIENT  =  .08687     PRESSURE  COEFFICIENT  =  .19671 
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APPENDIX  C 
CALIBRATION  OF  THE  FLOW  METER 
Since  flow  rates  were  very  small  without  obstructions  in  the 
discharge  line,  it  was  felt  that  some  method  of  determining  flow  rates 
with  a  minimum  obstruction  would  best  suit  the  purposes  of  this  experi- 
ment c  The  only  obstruction  introduced  into  the  flow  in  the  flow  meter 
devised  was  a  Kiel  probe.  For  purposes  of  calibration  the  discharge 
pipe  was  connected  to  an  air  supply  line  in  which  was  installed  a 
standard  sharp-edged  orifice  type  flow  meter,  so  that  all  the  air  which 
passed  through  the  standard  meter  also  passed  through  the  discharge  pipe0 
The  static  pressure,  Pg,  the  difference  between  the  centerline  total 
pressure  and  the  static  pressure, AP1 ,  and  the  centerline  total  tem- 
perature, T2,  were  measured  in  the  discharge  pipe.  Fig.  5.  The 
following  two  equations  were  to  be  used  for  determining  flow  rate. 

(1)  P^V2  a      fAP' 

2gRT2 

(2)  M  -  P^  AV 

RT2 

combining, 

(3)  w  =  A(2AP»fgP;/RT2)2 

Everything  can  be  measured  except  the  factor  T  .  The  object  of  the 

calibration  then  is  to  determine  T  .  This  factor  was  not  a  constant  so 

following  the  usual  practice  in  flow  meters,  "J  was  plotted  as  a  function 

of  Reynolds  number.  Fig.  18.  In  this  graph  a  modified  Reynolds 

number,  Rex,  was  used,  where, 

Rex  =  Re/3600  =  VD/» 

3600/1 
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APPENDIX  D 
APPLICATIONS  OF  THE  HOT  WIRE  ANEMOMETER 
The  hot  wire  anemometer  operates  on  the  principle  that  the  resist- 
ance of  most  metal  wire  is  a  linear  function  of  the  temperature ; 

R  =  %(J.  fo^d-  T0)] 
where,  RQ  is  the  resistance  at  some  reference  temperature,  TQ;  and  c<0 
is  the  -  temperature  coefficient  of  resistivity  at  that  temperature. 

It  has  also  been  found  that,  if  a  very  fine  wire  is  heated  by  an 
electric  current  and  placed  in  a  moving  gas  stream,  the  square  of  the 
current  required  to  maintain  a  constant  resistance  in  the  wire  is  a 
linear  function  of  square  root  of  the  product  of  the  static  pressure 
and  the  effective  cooling  velocity  of  the  gas  for  Mach  Nos.  less 
than  0.3. 

I2oc(pVe)* 
If  ft  is  the  angle  between  the  wire  axis  and  the  flow  direction  of 
the  cooling  gas,  then: 

Ve2  V  sin/*  A5°  ^  /3  *-   135° 

or     V  ^  V  sin£  (l  -  1.2(D/L)£  cos2)3  )    5°  <=  /S  <  175° 

.0002  ^  D/L  <   .001 
where,  D/L  is  the  ratio  of  diameter  to  length  of  the  wire  and  V  is  the 
magnitude  of  the  velocity  of  the  cooling  gas. 

Thus  if  an  air  stream  of  known  velocity  and  direction  is  available, 
a  calibration  curve  may  easily  be  obtained  by  taking  a  current  reading 
with  the  hot  wire  in  still  air  and  another  with  the  wire  in  the  air 
stream.  The  two  points  may  then  be  plotted,  connected  with  a  straight 
line  and  used  as  a  calibration  curve.  Fig.  23. 
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(Pve)i 

Fig.  23.  Sample  Calibration  Curve 
A  Wheatstone  bridge  is  normally  used  to  observe  the  resistance  in 
the  hot  wire  and  maintain  it  at  a  constant  value.  Fig.  24.. 


Fig.  24-.  Anemometer  Circuit 

Once  a  calibration  curve  is  available  the  hot  wire  may  be  used  to 
determine  velocities  of  airstreams.  It  will  be  noted  that  for  measuring 
velocity  with  the  hot  wire,  a  J3  of  90°  should  be  used.  Ve  is  then  taken 
as  the  velocity  because  in  this  orientation  the  measured  velocity  is 
very  insensitive  to  directional  orientation  of  the  wire. 

If  the  current  is  held  constant  and  the  hot  wire  immersed  in  the 
moving  gas  stream,  the  resistance  will  vary  as  a  function  of  the  product 
of  the  static  pressure  and  the  effective  cooling  velocity  of  the  gas. 
Fig.  25. 


A? 


pV 
r  e 

Fig.  25.  Sample  Characteristics  Curve 
This  fact  is  useful  in  determining  flow  direction.  For  example, 
at  a  current  of  1^,  as  jB  is  varied  through  360°,  the  resistance  varies 
as  indicated  by  the  heavy  lined  portion  of  the  above  curve  while  the 
wire  is  immersed  in  a  constant  velocity,  constant  static  pressure  stream 
of  gas.  The  minimum  value  of  R  would  occur  when^Or  90°,  270°  and  the 
maximum  R  would  occur  when  /3  _  0°,  ISO0.  Another  useful  observation 
is  that  the  resistance  would  be  the  same  value  for  /6  ■  10°  as  for 
J&  =  350°.  The  galvanometer  in  the  Wheatstone  bridge  will  indicate 
the  resistance  of  the  wire. 

It  will  be  noted  that  the  directional  sensitivity  of  the  hot  wire 
is  greatest  for   near  zero  for  two  reasons: 

(1)  IdR  I  is  greater  for  lower  Ve 
|dVel 

(2)  dVe   is  maximum  near  /Q  =  0 

w 

To  measure  the  flow  direction  no  calibration  is  required  except  to 
determine  the  directional  orientation  of  the  hot  wire  axis  with  respect 
to  some  device  for  measuring  angular  displacement  of  the  axis  such  as 
a  protractor. 

The  hot  wire  elements  are  normally  mounted  on  probes  across  two 
metal  posts  which  are  of  a  material  which  is  readily  soldered.  Fig.  26  a 
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To  Wheatstone 
Bridge 


Hot 

Wire 


Fig.  26.  Hot  Wire  Probe 

This  arrangement  fixes  the  wire  axis  approximately  perpendicular 
to  the  probe  axis.  A  protractor  may  be  fixed  to  the  probe  in  such  a 
manner  that,  as  the  probe  is  rotated  about  its  axis,  the  angular  dis- 
placement from  some  reference  position  may  be  measured.  By  suitable 
calibration,  any  given  orientation  of  the  wire  axis  may  be  made  to 
correspond  to  a  known  reading  on  the  protractor. 

Using  this  instrument  the  flow  direction  may  be  determined  in  the 
plane  in  which  the  wire  is  turned.  To  accomplish  this,  the  probe  is 
inserted  into  an  airstream  so  that  the  probe  axis  is  approximately  per- 
pendicular to  the  flow  direction  and  rotated  about  its  own  axis.  The 
resistance  of  the  hot  wire  will  vary,  the  maximum  value  being  when  the 
wire  axis  is  in  alignment  with  the  component  of  velocity  in  the  plane 
of  the  investigation.  As  an  example,  suppose  the  maximum  resistance 
is  observed  when  the  protractor  reads  approximately  34-5°.  It  is 
difficult  to  tell  exactly  when  the  maximum  resistance  is  encountered 
because  the  galvanometer  needle  is  fluctuating.  However,  when  the  pro- 
tractor reads  360°,  the  galvanometer  needle  fluctuates  about  a  mark  which 
will  indicate  some  slightly  smaller  resistance.  The  needle  fluctuates 
about  this  same  mark  when  the  protractor  reads  331.2°.  The  protractor 


u 


reading  for  a  maximum  resistance  should  then  be  halfway  between  the 
two  readings,  or  34-5 . 6°. 

It  may  be  noted  that  it  is  not  necessary  for  the  two  support  posts 
to  be  the  same  length,  Fig.  26,  nor  that  the  wire  axis  be  perpendicular 
to  the  probe  axis.  The  directional  sensitivity  suffers  when  the  wire 
axis  is  not  perpendicular  to  the  probe  axis,  however. 

One  possible  way  of  eliminating  the  effect  of  perturbation 
velocities  would  be  to  use  a  probe  with  two  hot  wires  mounted  in  a  cross 
fashion.  Fig.  27.  The  two  wires  could  be  placed  in  different  legs 
of  the  bridge  so  that  effects  of  perturbations  would  be  cancelling. 
After  calibration  in  a  known  direction  stream  it  could  be  used  for 
direction  determination.  This  instrument  would  also  have  double  the 
sensitivity  of  a  single  wire. 


Fig.  27.  Scheme  Fefi  Using  A  Two  Wire 
Probe  For  Measuring %low  Direction. 
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APPENDIX  E 
THE  COMPUTER  PROGRAM 
A  computer  program  was  used  to  reduce  the  data  and  is  included 
herein . 

The  following  symbols  were  used  in  the  program: 
Inputs 

run  number 

date 

month  indicator:.,  *-l  March,  -1  April 

position  no.  where  data  was  taken 


N 

M 

L 

NN 

P2 


PST 

PATM 

DLP 

T2 

BP 

ALFA1 
ALFA2 
TARE 

Outputs 
PHI  2 
PSI2 
Y 
ANGLE 


height  of  water  column  in  tube  to  which  pressure  tap 
at  outlet  flange  is  attached  -  cm.  of  ^0 


height  of  water  column  in  tube  to  which  static  pressure 
line  of  flow  meter  is  attached  -  cm.  of  HO 

height  of  water  column  in  tube  open  to  atmosphere 
cm.  of  H20 

difference  between  total  and  static  pressure  in  flow 
meter  -  in.  ^0 

total  temperature  in  air  at  flow  meter  -  °C 

barometer  reading  -^in.  Hg 

one  of  the  angles  read  at  each  depth  -  degrees 

the  other  angle  read  at  each  depth  --.degrees 

the  negative  of  flow  angle  read  on  the  protractor 
at  position  #2  on  the  centerline  of  the  inlet  duct 

flow  coefficient 

pressure  coefficient 

depth 

flow  direction 
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Other 

WI  mass  flow  rate  at  beginning  of  ran 

WE  mass  flow  rate  at  end  of  run 

PSI2I  pressure  coefficient  at  beginning  of  run 

PSI2E  pressure  coefficient  at  end  of  run 

PHI2I  flow  coefficient  at  beginning  of  run 

PHI2E  flow  coefficient  at  end  of  run 

PRI  pressure  ratio  at  beginning  of  run 

PRE  pressure  ratio  at  end  of  run 

W  mass  flow  rate 

PR  pressure  ratio 

ALFA  intermediate  step  in  averaging  process 

Other  in 
Subroutine 

P2A  pressure  at  discharge  flange  of  pump  -  PSF 

p 

A2  discharge  area  of  impeller  -  ft 

U2  impeller  peripheral  speed  -  ft/sec 

R  gas  constant 

PSTA  flow  meter  static  pressure  -  PSF 

TT2  same  temperature  as  72  -  °R 

BPA  barometric  pressure  -  PSF 

RHO  density  at  flow  meter  -  lbm/ft3 

ZETA  flow  rate  factor 

T2F  same  temperature  aa  T2  -  °F 

Z  viscosity  of  air  in  flow  meter 

VAV  average  velocity  in  flow  meter  -  ft/sec 

A  cross -sectional  area  of  discharge  pipe  -  ft2 

REX  modified  Reynolds  number  in  discharge  pipe 

WREF  referred  flow  rate  ~  ft2 
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For  each  run,  this  program  requires  seven  data  cards.  Cn  the  first 
card,  the  values  of  N,  M,  L  and  NN  are  entered  in  four  fixed  point  fields 
of  ten.  The  second  and  third  cards  are  six  floating  point  fields  of 
ten  in  which  are  entered  the  values  obtained  for  P2,  PST,  DLP,  T2,  PATM, 
and  BP„  The  second  card  contains  those  values  measured  at  the  beginning 
of  the  run  and  the  third  card,  those  values  measured  at  the  end  of  the 
rum  The  next  four  cards  contain  eight  floating  point  fields  of  ten. 
The  first  of  these  contains  the  values  of  ALFAl(l)  through  ALFAl(8). 
The  second  contains  the  values  of  ALFA1(9)  through  ALFA1(15)  and  TARE. 
The  last  two  cards  are  the  same  as  cards  four  and  five,  except 
substitute  ALFA2  for  ALFA1. 

The  program  first  enters  the  subroutine  twice  to  compute  flow 
coefficient  and  pressure  coefficient  from  the  data  on  cards  two  and 
three.  It  first  changes  all  pressures  to  PSFA  and  temperature  to  °R. 
If  different  manometers  and  thermometers  are  used,,  some  of  the  constants 
may  require  change .  It  then  does  two  iterations  to  find  ¥  because  of 
the  interdependence  of  ¥  and  ZETA.  The  formula  for  ZETA  in  terms  of 
REX  should  not  be  used  for  REX  greater  than  70 .  Fig.  17.  The 
subroutine  then  computes  the  pressure  coefficient  and  the  flow 
coefficient. 

The  flow  coefficients  and  pressure  coefficients  at  the  beginning 
and  end  of  the  run  are  then  averaged.  The  first  part  of  the  table  is 
then  printed.  Next  the  two  angles  measured  for  each  depth  are  averaged 
and  printed  with  its  corresponding  depth. 
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20 
25 


21 

22 
23 


PROGRAM 
DIMENSI 
Y( 1  )  =  0 
DO  20  J 
A  =  J  - 1 
Y(J)  = 
Y(  15)  = 
READ  99 
CALL  FL 
IF(PRI ) 
CALL  FL 
W=( WE+W 
PR=(PRL 
PSI2=(P 
PHI2=(P 
PRINT  1 
PRINT  1 
PRINT  1 
IF(L)21 


PREROT 
ON  Y(  15) 
,.  1 
=  2,14 


,  ALFAK  15),  ALFA2(  15)  ,  ANGLE (  15) 


PRINT 
GO  TO 
PRINT 
PRINT 
PRINT 
PRINT 
READ 


1 
2 
1 
1 
1 
1 
30 


0.5*  »  A 

6.9 
,  N,M,L, 
OW  ( Wit P 
7,7,6 
OW  (WE,PSI2.E,P 
I  )/2.0 
+PRI )/2. 
SI2I+PSI 
HI2I+PHI 
10 

11,  N 
12 

,21,22 
21.i'N«M 
3 

22,N,M 
NN 
PHI2,PSI2 


NN 
SI2I,P 


0 

2E)/2. 

2E)/2. 


HI2I,PRI  ) 
HI2E,PRE) 


1 
2 

3 

5 

k 

30 

110 

111 

112 

121 

122 

123 

1240 

1 

125 

200 

99 

300 

7 


READ  30 

DO  30  K 

ALFA=  A 

IF(ALFA 

IF(ALFA 

ALFA=AL 

GO  TO  4 

IF( ALFA 
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